o 

(N 



p 

o 



Quark-hadron phase transition in Brans-Dicke brane 

gravity 

K. Atazadeh^*, A. M. Ghezelbash^t and H. R. Sepangi^^ 
^Department of Physics, Shahid Beheshti University G. C, Evin, Tehran 19839, Iran 
^Department of Physics and Engineering Physics, University of Saskatchewan, 
Saskatoon, Saskatchewan S7N 5E2, Canada 

January 13, 2013 



Abstract 



A standard picture in cosmology has been emerging over the past decade in which a phase 
' ^ ' , transition, associated with chiral symmetry breaking after the electroweak transition, has occurred 

at approximately 10~^ seconds after the Big Bang to convert a plasma of free quarks and gluons 
^ , into hadrons. In this paper, we consider the quark-hadron phase transition in a Brans-Dicke 

bJ)' brane world scenario within an effective model of QCD. We study the evolution of the physical 

quantities relevant to quantitative description of the early universe, namely, the energy density, 
temperature and the scale factor before, during, and after the phase transition. We show that for 
^ , different values of the Brans-Dicke coupling, w, phase transition occurs and results in decreasing 

CO ' the effective temperature of the quark-gluon plasma and of the hadronic fluid. We then move on to 

, consider the quark-hadron transition in the smooth crossover regime at high and low temperatures 

' and show that such a transition occurs and results in decreasing the effective temperature of the 



quark-gluon plasma during the process of quark-hadron phase transition. 

1 Introduction 



Over the past decade the possibility that the observable universe is a brane [1] embedded in a higher 
^ , dimensional space-time has been the hallmark of great many research works. This scenario has 

motivated intense efforts to understand the case where the bulk is a 5-dimensional anti de-Sitter 
space. In this setup, gravitons are allowed to penetrate into the bulk but are localized on and around 
the brane [2] . It was then shown that in a background of a non-factorizable geometry an exponential 
warp factor emerges which multiplies the Poincare invariant 3-1-1 dimensions in the metric. The 
existence of branes and the requirement that matter fields should be localized on the brane lead to a 
non-conventional cosmology which has seeded a large number of studies. Of interest in the present 
study are the brane-world models in the context of Brans-Dicke (BD) gravity. Interestingly, it has 
been shown that in such a scenario and in the presence of a BD field in the bulk [3] the conservation 
equation for the matter field can be satisfied. It would therefore be of interest to study the quark- 
hadron phase transition in the context of BD brane world theory. The question of quark-hadron 
phase transition in the context of conventional brane-world models has been addressed previously 
[4, 5]. 

Standard cosmology suggests that as the early universe expanded and cooled, it underwent a 
series of symmetry-breaking phase transitions, causing topological defects to form. It is the study of 
such phase transitions that would pave the way for a better understanding of the evolution of the 
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early universe, characterized by the existence of a quark-gluon plasma undergoing a phase transition. 
In what follows we focus attention on possible scenarios which might have occurred to allow such a 
phase transition come to the fore. We generally follow the discussion presented in [4] which puts the 
quark-gluon phase transition in a cosmologically transparent perspective. 

The existence of phase transition from the quark-gluon plasma phase to hadron gas phase is a 
definite prediction of QCD. However, the phase transition in QCD can be characterized by a truly 
singular behavior of the partition function leading to a first or second order phase transition, and 
also it can be only a crossover with rapid changes in some observables, strongly depending on the 
values of the quark masses. The possibility of a phase transition in the gas of quark-gluon bags was 
demonstrated for the first time in [6]. Most studies have shown that one may obtain first, second 
and higher order transitions. In addition, the possibility of no phase transitions was pointed out in 
[7]. Recently, lattice QCD calculations performed for two quark flavors suggest that QCD makes a 
smooth crossover transition at a temperature of Tc ~ 150 MeV [8]. Such a phase transition could 
be responsible for the formation of relic quark-gluon objects in the early universe which may have 
survived. In this paper, our study of phase transition is based on the ideas proposed in the first 
reference in [6], where it was shown that under certain conditions a gas of extended hadrons could 
produce phase transitions of the first or second order, and also a smooth crossover transition that 
might be qualitatively similar to that of lattice QCD. 

The cooling down of the color deconfined quark-gluon plasma below the critical temperature 
believed to be around ~ 150 McV makes it energetically favorable to form color confined hadrons 
(mainly pions and a tiny amount of neutrons and protons, since the net baryon number should be 
conserved). However, such a new phase does not form promptly. Generally speaking, a first order 
phase transition needs some supercooling to achieve the energy used in forming the surface of the 
bubble and the new hadron phase. A brief account of a first order quark-hadron phase transition in the 
expanding universe may be envisioned as follows [9]. When a hadron bubble is nucleated, latent heat is 
released and a spherical shock wave expands into the surrounding supercooled quark-gluon plasma. 
The plasma thus formed is reheated and approaches the critical temperature, preventing further 
nucleation in a region passed by one or more shock fronts. Bubble growth is generally described by 
deflagrations where a shock front precedes the actual transition front. The stopping of the nucleation 
occurs when the whole universe has reheated T^. The prompt ending of this phase transition, in 
about 0.05 /iS, renders the cosmic expansion completely negligible over this period. Afterwards, 
the hadron bubbles grow at the expense of the quark phase and eventually percolate or coalesce. 
Eventually, when all quark-gluon plasma has been converted into hadrons, neglecting possible quark 
nugget production, the transition ends. The physics of the quark-hadron phase transition and its 
cosmological implications have been extensively discussed in the framework of general relativistic 
cosmology in [10]-[23]. 

As is well known, the Friedmann equations in branc-world scenarios differs from that of the stan- 
dard AD cosmology in that they result in an increased expansion rate at early times. We expect this 
deviation from the standard AD cosmology to have noticeable effects on the cosmological evolution, 
especially on cosmological phase transitions. In the context of brane-world models, the first order 
phase transitions have been studied in [24] where it has been shown that due to the effects coming 
from higher dimensions, a phase transition requires a higher nucleation rate to complete and baryoge- 
nesis and particle abundances could be suppressed. Recently, the quark-hadron phase transition was 
studied in a Randall-Sundrum brane-world scenario [4]. Within the framework of first order phase 
transitions, the authors studied the evolution of the relevant cosmological parameters (energy density, 
temperature, scale factor, etc.) of the quark-gluon and hadron phases and the phase transition itself. 
In another attempt, the phase transition of quarks and gluons was studied in a brane-world model 
in which the confinement of matter fields on the brane is achieved through a confining potential [5]. 
As was mentioned above, it would therefore be of interest to study phase transitions of this nature 
in the context of a BD brane-world scenario and this is what we intend to do in what follows. In 
addition, as recent calculations in lattice QCD strongly favor a smooth crossover transition [8], the 
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study of the formation of hadrons in such a scenario seems to be of particular importance for a better 
understanding of the subtleties of the early universe and this is what we shall present in the last 
section. 



2 Field equations in Brans-Dicke brane scenario 

We start by writing the action for the Brans-Dicke brane-world [25] 



(1) 



where TZ is the Ricci scalar associated with the 5-dimcnsional space-time metric g^^g, is a scalar 
field which we shall call the BD field, w is a dimensionless coupling constant which determines the 
coupling between gravity and the BD scalar field and >Cm represents the Lagrangian for the matter 
fields. Latin indices denote 5-dimensional components (^4, B = 0, . . . , 5) and for convenience we 
choose = 87rG(5) = 1. The variation of the action with respect to g^g and 4> yields the field 
equations 



Gab = T^AB — -^Bab'^ 



1 



Tab + Tab 



where 
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+ [<1>;AB-9aB<1>' -,0] , 



(2) 



(3) 



and T = T^c is the trace of the energy-momentum tensor of the matter content of the 5-dimensional 
space-time. Note the factor 3a; -I- 4 in the denominator on the right hand side of the BD field equation 
instead of the familiar 2c<j + 3 in the 4-dimensional case [26]. This is determined by requiring the 
validity of the equivalence principle in our setup, see [27] for a discussion of this topic in the context 
of 4-dimensional BD theory. 

Being interested in cosmological solutions, we take the spatially flat cosmology {k = 0) and 
consider a 5-dimensional fiat metric of the following form 

ds'^ = -n^{T, y)dT^ + a^{T, y)5ijdx^dx^ + 6^(r, y)dy^ , (5) 

where i,j = 1,2,3. We also assume an orbifold symmetry along the fifth direction y — )• —y. Next we 
define the energy-momentum tensor 

^^^B = T'^Blbulk + T'^Bjbrane , (6) 

where the subscripts "brane" and "bulk" refer to the corresponding energy-momentum tensors. For 
simplicity we assume that the bulk is devoid of matter other than the BD scalar field. The brane 
matter field is held at y = with the following energy momentum tensors 

r^B|brane = ^diag(-/9,p,p,p,0) and T^B [bulk = diag(0, 0, 0, 0, 0), (7) 

where p = + X and p = p'' — X. The above expressions are written assuming that the brane has 
ordinary matter with tension A and that the bulk is empty. There are several constraints suggested 
for the brane tension A. One is that suggested by the big bang nucleosynthesis, A > 1 MeV^ [28]. A 
much stronger bound for A is due to the null results of submillimeter tests of Newton's law, giving 
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A > 10^ GeV^ [29]. An astrophysical lower limit on A which is independent of the Newton's law and 
cosmological limits has been studied in [28] , leading to the value A > 5 x 10^ MeV^ which is constraint 
we will be using in our model. 

Using metric (5) we are now able to write the equations of motion. The (0, 0) component reads 



a I a b 
a \ a b 



V? fa" 
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a' fa' 



+ 
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where 
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The (0, 5) component takes the form 
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Finally, for the (5, 5) component one has 
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The equation of motion for the BD field reads 



■ I „d b h 
+ </. 3- + --- 
\ a b n 



a' b' n' 
3--- + - 
a b n 



-n 



3a; + 4' 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



where a dot represents the time derivative with respect to r and the prime corresponds to derivatives 
with respect to y. Note that in the above derivation we have assumed (p = 4>ij-,y)- We make 
the assumption that the metric and the BD field are continuous across the brane localized at y = 0. 
However, their derivatives can be discontinuous at the brane position in the y direction. This suggests 
the second derivatives of the scale factor and the BD field will have a Dirac delta function associated 
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with the positions of the brane. Since the matter is localized on the brane it will introduce a delta 
function in the Einstein equations which will be matched by the distributional part of the second 
derivatives of the scale factor and BD field. For instance at y = 0, we have [25] 



/" = /" + [/']%), 



(17) 



where the hat marks the non-distributional part of the second derivative of the quantity. The part 
associated with a delta function, [/'], is a jump in the derivative of /. Here / could be any of the 
three quantities a, n or 0. The jump in / at y = can be written as 



[/'] = /'(o+)-/'(o-), 

and the mean value of the function / at y = is defined by 



m 



/'(o+)-/'(o-) 



(18) 



(19) 



After substituting equation (17) in the Einstein field equations it is possible to find the jump 
conditions for a and n by matching the Dirac delta functions appearing on the left-hand side of the 
Einstein equations to the ones coming from the energy-momentum tensor, equation (6). For the BD 
field one has to use equation (16) to evaluate the jump conditions. We therefore find 



n 



where 



(3w + 4)(/.„ 
1 



(3a; + ^)(t>^ 
2 



3(w + l)p + (2a; + 2,)p 



IP, 



(20) 
(21) 
(22) 



7 = l^i'^Wm - 1) 



(23) 



with Wm = ^ and the subscript stands for the brane at y = 0. The first two conditions, equations (20) 
and (21), are equivalent to Israel's junction conditions in general relativity (see [30] for a discussion 
of its application in the context of brane- wor Ids) . It is important to note that the above jump 
conditions at y = depend on the energy density and pressure component of the brane and the 
induced curvature on the brane. Interestingly, for the radiation dominated phase on the brane, 
p = 3p, the jump condition for cp does not vanish and is proportional to the energy density and 
pressure component of the induced curvature on the brane. 

Using the (0, 0) component of the Einstein field equations for the brane located at y = and the 
equations representing the jump conditions (20), (21) and (22) one gets the Friedmann equation as 
follows 



4(30.' + 4)2(^2 



I (3p - p)2 + (2 + 3a; + uj^)p'^ - ojpp - 2p'^ 



(24) 



where H 



p^ + X and p — 



p" 



A, 



and T = Note that p and p consist of two parts, that is p ■ 
where A is tension on the brane. Using the (0, 5) component of the Einstein equation and substituting 
equations (20) and (21) we get the continuity equation for the matter on the brane 



p + 3{p + p)H = 0. 



(25) 



While deriving the above equations we have assumed that, from the point of view of the brane 
observer, the extra dimension is static, that is 6 = 6o. We have also fixed the time in such a way that 
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no = 1, corresponding to the usual choice of time in conventional cosmology. Taking the mean value 
of the BD field equation we obtain an equation of motion for (j) on the brane 

^ + SHT = ^^P^. (26) 

(/.Q (3w + 4)2,/.2 

Note that in order to obtain equation (26) wc also have to assume that the non-distributional part 
of 4>" vanishes, otherwise, a term involving 0" will appear in the BD field equation. As we shall 
see in the next section it is possible to obtain cosmologically interesting solutions which verify this 
condition. Equations (24) and (25) are dynamical equations in our BD brane scenario describing the 
evolution of the universe. In the next section we shall examine these equations for the quark-hadron 
phase transition in the early universe. 



3 Quark-hadron phase transition 



The quark-hadron phase transition is a notion fundamental to the study of particle physics, partic- 
ularly in the context of lattice gauge theories. However, it is an integral part of any study dealing 
with the underlying mechanisms responsible for the evolving universe at its early stages of formation 
in which a soup of quarks and gluons interact and undergo a phase transition to form hadrons. It is 
therefore essential to have a overview of the basic ideas before attempting to use the results obtained 
from such a phase transition and apply them to the study of the evolution of the early universe within 
the context of the BD brane-world scenario. In this regard, a well written and concise review can be 
foTind in [4] and the interested reader should consult it. Here, it would suffice to mention the results 
relevant to our study and leave the details of the discussion to the said reference. 

We start from the equation of state of matter in the quark phase which can generally be given in 
the form 



pl=3a^T^ + V{T), 



p'=a^T'-V{T), 



(27) 



where = {■K'^/90)g^, with = 16 + {21/2)Nf + 14.25 = 51.25 and Np = 2 with V{T) being the 
self-interaction potential. For V{T) we adopt the expression [21] 



(28) 



where B is the bag pressure constant, ax = 77r2/20 and 7^ = ^2/4 with the mass of the strange 
quark in the range nig G (60 — 200) McV. In the case where the temperature effects are ignorablc, 
the equation of state in the quark phase takes the form of the MIT bag model equation of state, 
^9 = (p^ — 4i?)/3. Results obtained in low energy hadron spectroscopy, heavy ion collisions and 
phenomenological fits of light hadron properties give 

5^4 between 100 and 200 MeV [31]. 
Once the hadron phase is reached, one takes the cosmological fluid with energy density pj^ and 
pressure p'^^ as an ideal gas of massless pions and nucleons obeying the Maxwell-Boltzmann statistics. 
The equation of state can be approximated by 



1 



(29) 



where a,r = {'^'^/^0)9h and = 
P^Tc) = p^iTc) [9], and is given by 



17.25. The critical temperature Tc is defined by the condition 



1T + ^Jit + 4-B(a^ +aT- a^) 
2{a +aT- a^) 



1/2 



(30) 



If we take nis = 200 MeV and S^/^ = 200 MeV, the transition temperature is of the order Tc ~ 125 
MeV. It is worth mentioning that since the phase transition is assumed to be of first order, all the 
physical quantities exhibit discontinuities across the critical curve. 
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4 Behavior of Brans-Dicke brane universe during quark-hadron phase 
transition 



We are now in a position to study the phase transition described above. The framework wc are 
working in is defined by the BD brane-world scenario for which the basic equations were derived in 
section 2. The physical quantities of interest through the quark-hadron phase transition are the energy 
density p, temperature T and scale factor Og. These parameters are determined by the Priedmann 
equation (24), conservation equation (25) and the equations of state, namely (27), (28) and (29). To 
start, we consider the evolution of the BD brane-world before, during and after the phase transition 
era. 



4.1 Behavior of temperature 

Let us consider the era preceding the phase transition for which T > Tc and the universe is in the 
quark phase. Use of equations of state of the quark matter and the conservation of matter on the 
brane, equation (25), leads to 

ao 3aq T 6aqT^' ^ ' 

Integrating the above equation immediately gives 

aJT)^cT exp(-^J^), (32) 

where c is a constant of integration. 

To consider the phase transition in the BD brane model, for simplicity, we take an ansatz in the 
form of a relation between scale factor on the brane (a,,) and the BD scalar field (^o) as follows 

</.„(r) = K(^)> (33) 
where n and n are constant. Thus, using equation (33), the Priedmann equation (24) can be written 



as 



I (3p - pf + (2 + 3w + uj^)p^ - ppco - 2p2] . (34) 



2[(3a; + 4)'(2 + 2n-|n2)](/)2 

We may now proceed to obtain an expression describing the evolution of temperature of the BD 
brane universe in the quark phase by combining equations (27), (28), (31), (32), (33) and (34), leading 
to 



dr /x(3a; + 4){A^T^ + A^)^2{2 + 2n- fn^) ^"""^ V 



1/2 

(35) 



^(ViT) - Xf + (2 + 3a; + c^^)(pj(r) + A)^ - ^^(^^(r) - X)ipl{T) + A) - 2(pJ(r) - A)^ 



where we have denoted 

Ao = l-^, (36) 
6a, 

Equation (35) may be solved numerically and the result is presented in figure 1 which shows the 
behavior of temperature as a function of the cosmic time r in a BD brane-world filled with quark 
matter for different values of u with n = 0.05 , fx = 2 x 10^ and A = 10 x 10^ MeV^. 
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Figure 1: The behavior of T(r) as a function of time (t) for /.i = 2 x lO'^, A = 10 x 10* MeV*, n = 0.05 and different 
values of a;: a; = 1 x 10^ (solid curve), ui = 2.3 x 10^ (dashed curve), uj = 2.4 x 10^ (dotted curve) and co = 2.5 x 10^ 
(dotted-dashed curve). We have taken B^/** = 200 MeV. 



4.2 Temperature behavior with V(T) = B 

One may gain considerable insight in the evolution of cosmological quark matter in our BD brane- 
world by taking the simple case in which temperature corrections can be neglected in the self inter- 
acting potential V . Then V = B = const, and equation of state of the quark matter is given by that 
of the bag model, namely = — AB^ /3. Equation (25) may then be integrated to give the scale 
factor on the brane as a function of temperature 

«„(r) = |, <i,,=^Jic^T-\ (37) 

where c is a constant of integration. 

Using equations (34) and (37), the time dependence of temperature can be obtained from equation 

dT -r"+' 



dr /i(3a; + 4)y/2(2 + 2 



n — Trn 




-f{B - Xf + {2 + 3u: + u^){p\{T) + A)2 - u; '-^^ a (pj(r) + A) - 2 '-^^-^ 



where we have taken c = 1 and (T) is given by 

p\{T) = 2,aJ^ + B. (39) 

The above equation can be solved numerically and the result is given in figure 2 which shows the 
behavior of temperature as a function of cosmic time r in a BD brane- world filled with quark matter 
for different values of uj with p = 2 x 10^ , n = 0.05 and A = 10 x 10^ MeV^ as the self-interacting 
potential, V{T), is considered to be a constant. 

4.3 Formation of hadrons 

During the phase transition, the temperature and pressure are constant and quantities like the entropy 
S = sa^ and enthalpy W = [p + p) are conserved. Also, during the phase transition /o''(i) decreases 
from Pg{Tc) = pq to p\{Tc) = p^ ■ For phase transition temperature of Tc = 125 MeV we have 
Pq ^ 5 X lO^MeV* and pn ~ 1-38 x 10^ MeV^, respectively. For the same value of the temperature 
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Figure 2: The behavior of T(t) as a function of time (r) for = 2 x lO'^, A = 10 x 10** MeV*, n = 0.05 and different 
values of a;: ui = 100 (solid curve), a; = 1 x 10'^ (dashed curve), lo = 2 x 10^ (dotted curve) and u> = 2.2 x 10^ 
(dotted-dashed curve). We have taken B^^'^ = 200 MeV. 



the value of the pressure of the cosmological fluid during the phase transition is ~ 4.6 x 10^ MeV^. 
Following [9, 4], we replace (t) by /i(t), the volume fraction of matter in the hadron phase, by 
defining 

(r) = p^h{T) +Pq[1- h{T)] =Pq[1 + mh{T)] , (40) 

where m = (p^ ~ Pq) / Pq- The beginning of the phase transition is characterized by /i(rc) = where 
Tc is the time representing it and (tc) = Pq , while the end of the transition is characterized by 
h{Th) = 1 with Th being the time signaling the end and corresponding to p'' (r/i) = p^. For t > Th 
the universe enters into the hadronic phase. 
Equation (25) now gives 

«o _ 1 {Ph-Pq)^ _ 1 rh ^^^^ 
«o ^Pq+Pc + {Ph - Pq)^^ 3 1 + r/i' 

where we have denoted r = (p^ ~ Pq) / {Pq + Pc) ■ The relation between the scale factor on the brane 
and the hadronic fraction /i(r) may now be obtained from the above equation 

a„(r) = a, (r^) [1 + rh{T)]-'/' , = ^ (^e) [1 + rhir)]-"/' , (42) 

where use has been made of the initial condition h (tc) = 0. Now, using equations (34) and (42) we 
obtain the time evolution of the matter fraction in the hadronic phase 

dh 3[1 + rh{t)]^ 1 



dr pr{3uj + 4)aJ^ (tc) ^^2(2 + 2n - | 



^(^P. - X{t)p^ - 4A)2 + (2 + 3u; + u^){x{t)p^ + Xf - u{p^ - A)(x(t)Pg + A) - 2(p, - A)^ 



(43) 

Here we take a"(Tc) = 1 and x(0 = 1 + rnh{t). Figure 3 shows variation of the hadron fraction /i(t) 
as a function of r for different values of oj with p = 20000 , n = 0.05 and A = 8 x 10* MeV^. 
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Figure 3: The behavior of ft(r) as a function of time (r) for fi = 20000, A = 8 x 10* MeV, n = 0.05, V{T) = B 
and different values of tj: cj = 1 x 10^ (soUd curve), oj = 2.1 x 10^ (dashed curve), uj — 2.3 x 10^ (dotted curve) and 
uj = 2.47 x 10^ (dotted-dashed curve). We have taken B^/* = 200 MeV. 



T(T) 
130 r 




Figure 4: The behavior of T{t) as a function of time (r) for /i = 20000, A = 5 x 10* MeV, n = 0.02 and different 
values of oj: a; = 15 x 10^ (solid curve), a; = 13 x 10^ (dashed curve), ui = 9 x 10^ (dotted curve) and cj = 1 x 10^ 
(dotted-dashed curve). We have taken B^^* = 200 MeV and Tc = 125 MeV. 



4.4 Pure hadronic era 

After the phase transition, the energy density of the pure hadronic matter is given by = 3p^ = 
SottT^. The conservation equation on the brane (25) leads to 

a,{T) = a,{th)TjT, 4>o = Ki^hWc/TT . (44) 

The temperature dependence of the Brans-Dicke brane universe in the hadronic phase is governed 
by the equation 

dT ^ r"+i [f^A^ + {uj^ + 1^ + f )(a,T4 + A)^] ^ 
dr ~ n[Tc%{th)Y{'ioj + 4)^/2(2 + 2n - |n^) 

Variation of temperature of the hadronic fluid filled BD brane universe as a function of r for different 
values of a; with /x = 20000 , n = 0.02 and A = 5 x 10^ MeV^ is represented in figure 4. 
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Figure 5: Left, the behavior of T(r) as a function of time (t) for p = 2 x lO'^, A = 10 x 10® MeV*, n = 0.05 and small 
and large values of a;, respectively : tj = 1 x 10° (solid curve) and uj = 2.518 x 10^ (dotted-dashed curve). Right, the 
behavior of /i(r) as a function of time (r) for = 20000, A = 8 x 10* MeV, n = 0.05, V'(r) = B and small and large 
values of uj, respectively: cj = 1 x 10^ (sohd curve) and oj = 2.5 x lO'' (dotted-dashed curve). We have taken B^^"^ = 200 



4.5 Effects of the BD coupling on phase transition 

It is important to emphasize that the BD coupling, w, plays an essential role in the model at hand 
and it is thus appropriate at this point to consider the effects of uj during the phase transition. 

We know that in the early universe the energy density is extremely high and thus in the high 
density regime the Hubble function in the DB brane world is proportional to both the energy density 
of the cosmological matter and oj, as can be seen from terms in the square bracket in equation (24). 
Moreover, in the early universe the BD coupling, which is an indication of the strength of the DB 
scalar field, is expected to play an important role and therefore its numerical value would characterize 
its influence. In general, in the limit w — )• oo, we recover the FRW equation in the Randall-Sundrum 
(RS) model [30] from equation (24). 

In the context of the BD brane world model, we have found that the temperature evolution of the 
universe is different from that of the RS brane world model. The temperature of the early universe 
in the quark phase is higher in the BD brane world scenario, as can be seen from the left graph in 
figure 5, where the dotted-dashed curve, to a high degree, corresponds to the RS brane world limit of 
large co. Hence a large value of the BD coupling, w, would significantly reduce the temperature of the 
quark-gluon plasma, and accelerates the phase transition to the hadronic era. Once the quark-hadron 
phase transition starts, the hadron fraction h is again strongly dependent on the BD coupling cj. From 
the right graph in figure 5 it is seen that for large values of uj, h[T) is much higher (dotted-dashed 
curve) than in the RS model and standard general relativity. The effect of an increase in oj on the 
brane is to strongly accelerate the formation of the hadronic phase and shorten the time interval 
necessary for the transition. A large uj tends to reduce the temperature of the hadronic fluid. From 
these figures it can be seen that for small values of cj, the rate of the phase transition is slower than 
the large values of uj. 

5 Lattice QCD phase transition 

As was mentioned in the introduction, lattice QCD calculations for two quark flavors suggest that 
QCD makes a smooth crossover transition at a temperature of ^ 200 MeV [8]. It is therefore 
necessary to have a brief review of the basic notions of this subject before using the results to study 
the universe at early times within the context of Brans-Dicke brane scenario without the cosmological 
constant. A, on the brane. 

Lattice QCD is an approach which allows one to systematically study the non-perturbative regime 
of the QCD equation of state. This approach has enabled the calculation of the QCD equation of 



MeV . 
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state using supercomputers [32] with two light quarks and a heavier strange quark on a {N^ = 6) x 32^ 
size lattice. The quark masses have been chosen to be close to their physical value, i.e. the pion mass 
is about 220 MeV . The equation of state was calculated at a temporal extent of the lattice Nt = 6 for 
which sizable lattice cut-off effects are still present [33]. The data for energy density p{T), pressure 
p{T) and trace anomaly p — 3p and entropy s, used in what follows are taken from reference [32]. We 
also note that besides the strange quark, one can also include the effect of the charm quark as well 
as photons and leptons on the equation of state. These have important cosmological contributions as 
was shown in [34]. Recent references on lattice QCD at high temperature can be found in [35]. 

Over the high temperature regime, radiation like behavior is seen as expected. However, in the 
region at and below the critical temperature Tc (~ 200 MeV ) of the deconfinement transition, the 
behavior changes drastically. This behavior change is also relevant for cosmological observables as we 
will see in the following. For high temperature, that is between 2.82 (lOOMeV) and 7.19 (lOOMeV), 
one can fit the data to a simple equation of state of the form 

piT) « a^^ (46) 
p{T) ^ aT"^. 

The values of a = 14.9702 ± 009997 and a = 4.99115 ± 004474 are found using a least squares fit [32]. 

While for times before the phase transition the lattice data matches the radiation behavior very 
well, for times corresponding to temperatures above Tc the behavior of the lattice data changes 
towards matter dominated behavior. We remark that lattice studies show that the QCD phase 
transition at its physical values is actually a crossover transition. 



5.1 High Temperature Regime 

Let us first consider the era before phase transition at high temperature where the universe is in the 
quark phase. Using the conservation equation of matter together with the equation of state of quark 
matter (46), one finds the following relation for the Hubble parameter 

JJ = 5 = -,^|, (47) 
a 3(a + a)T 

whose solution is given by 

a{T) = cT^(^, (48) 

where c is a constant of integration. 

One may now proceed to obtain an expression describing the behavior of temperature of the 
Brans-Dicke brane universe with respect to time in the quark phase. Using equations (34), (46), (47) 
and (48) one finds a differential equation for the temperature 



1 /2 

dT _ 3{a + a)T ( |(3p - p) + (2 + 3a; + uj^)p^ - upp - 2p'^ \ 

^" 4^0(30; + 4)r^ V 2[2(n + 1) - ^] ) 



(49) 



where we have set the constant c to unity. The transition region in the crossover regime can be 
defined as the temperature interval 282 MeV < T < 719 MeV. In view of the crossover nature of the 
finite temperature QCD transition, such definition is equivocal [36]. 

Equation (49) can be solved numerically and the result is plotted in figure 5 which shows the 
behavior of temperature of the Universe in the quark phase as a function of the cosmic time in 
Brans-Dicke brane cosmology for lj = 3 x 10^, in the interval 282 MeV < T < 719 MeV in the high 
temperature regime. We see that as the time evolves the universe becomes cooler. 
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Figure 6: The behavior of T{t) in the interval 282 MeV < T < 719 MeV as a function of r for /i 
and a; = 3 X lO'^ . 



10", n = 0.00004 



5.2 Low temperature regime 

Besides lattice QCD there are other approaches to the low temperature equation of state. In the 
framework of the Hadronic Resonance Gas model (HRG), QCD in the confinement phase is treated 
as a non-interacting gas of fermions and bosons [37]. The fermions and bosons in this model are 
the hadronic resonances of QCD, namely mesons and baryons. The idea of the HRG model is to 
implicitly account for the strong interaction in the confinement phase by looking at the hadronic 
resonances only since these are basically the relevant degrees of freedom in that phase. The HRG 
model is expected to give a good description of thermodynamic quantities in the transition region 
from high to low temperature [38]. The HRG result for the trace anomaly can also be parameterized 
by the simple form [36] 



e(r) _p-3p 



rp4^ 



rp4 



aiT + a2T^ + a^T + 



10 



(50) 
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with ai = 4.654 GeV-\ 02 = -879 GeV-^, 03 = 8081 GeV-"^, 04 = -7039000 GeV" 

In lattice QCD the calculation of the pressure, energy density and entropy density usually proceeds 
through the calculation of the trace anomaly Q{T) = p{T) — 3p{T). Using the thermodynamic 
identities, the pressure difference at temperatures T and Tiow can be expressed as the integral of the 
trace anomaly 



rp4 



low 



e(r'). 



(51) 



By choosing the lower integration limit sufficiently small, p{T\o„) can be neglected due to the ex- 
ponential suppression. Then the energy density p{T) = 6(T) -|- 3p(T) and the entropy density 
s{T) = [p + p)/T can be calculated. This procedure is known as the integral method [39]. Using 
equations (50) and (51) we obtain 



PiT) = r^T' + a,T' + ^T^ + ^T^ + ^T'\ 



(52) 



where 77 = —0.112. The trace anomaly plays a central role in lattice determination of the equation of 
state. The equation of state is obtained by integrating the parameterizations given in equations (50) 
over the temperature as shown in equation (51). 

Let us now consider the era before phase transition at low temperature where the universe is in 
the confinement phase and is treated as a non-interacting gas of fermions and bosons [37]. Using the 
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2.x 10 



4.X10-' 



6. X 10"' 



Figure 7: The behavior of T{t) in the interval 70 MeV < T < 180 MeV as a function of t ioY fi = 10^ 
and a; = 2 X lO". 



n = 0.00004, 



conservation equation of matter together with equation of state (52), one gets the following relation 
for the Hubble parameter 

^ d UrjT^ + 20aiT^ + A{T) . 



where 



a 3[4r/r4 + 5air5 + S(r)]' 



91 

AiT) = Ua2T^ + MasT^ + — T^^ 

5 

B{T) = la2T^ + 2a3T^ + la^T^*. 
6 5 



(53) 



(54) 



One can solve for the scale factor as 



a(r) 



T{75aiT + 35a2r3 + 30a3r4 + 2lT^o + 60??)i/3 ' 



(55) 



where c is a constant of integration. 

We can obtain an expression describing the behavior of temperature of the Brans-Dicke brane 
universe with respect to time in the quark phase. Upon using equations (34), (52), (53) and (55) one 
finds a differential equation for the temperature as follows 



dT 
d7 



3(4r?r4 + 5aiT^ + B) 



f (3p - p) + (2 + 3w + ^2)^2 _ _ 



1/2 



n{3oj + 4)a"(r)(12r/T3 + 20aiT^ + A) 



2[2(n + l) 



2 J 



(56) 

where we have set the constant c to unity. 

We can solve equation (56) numerically and the result is plotted in figure 6 which shows the 
behavior of temperature of the universe in the quark phase as a function of cosmic time in Brans- 
Dicke brane cosmology for uj = 2x 10^, in the interval 70 MeV < T < 180 MeV in the low temperature 
regime. 

The behavior of temperature with time can now be compared in the two regimes considered above, 
namely the first order phase transition and the smooth crossover transition. Figures 2 and 4 show 
the variation of temperature in the former while figurers 5 and 6 show such a variation in the latter. 
As can be seen, the general behavior is similar, that is, in both regimes the temperature drops as the 
time passes. We also note that the variation of T with cosmic time is somewhat different in the two 
approaches when considered in detail. In the smooth crossover regime where lattice QCD is used to 
investigate the high temperature behavior, we see that the slope is smooth relative to the first order 
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phase transition while at lower temperatures where HRG is used, the slope is steep compared to first 
order phase transition. Taking into account the energy range in which the calculations are done, one 
might conclude that these two approaches to the quark-hadron transition in the early universe do not 
predict fundamentally different ways of the evolution of the early universe. 

6 Conclusions 

In this paper, we have discussed the quark-hadron phase transition in a BD brane cosmological 

setting in which our universe is a thrcc-brane embedded in a 5-dimcnsional bulk space-time within an 
effective model of QCD. We studied the evolution of the physical quantities relevant to the physical 
description of the early universe; the energy density, temperature and scale factor, before, during, and 
after the phase transition. We found that for different values of oj phase transition occurs and results 
in decreasing the effective temperature of the quark-gluon plasma and of the hadronic fluid. We then 
compared our results with the results presented in [4] and [5]. In [4] the authors studied quark- 
hadron phase transition in a Randall-Sundrum brane model and showed that for different values of 
the barne tension A, phase transition occurs. Also in [5], the authors investigated the quark-hadron 
phase transition in a brane-world scenario where the localization of matter on the brane is achieved 
throTigh the action of a confining potential and showed that for different values of the parameters in 
their model, phase transition takes place. What has been lacking is a study of such phase transitions 
in the context of a BD brane world model, presented in this work, emphasizing the evolution of the 
physical quantities relevant to the physical description of the early universe. 

Finally, in section 5 we have considered the quark-hadron transition in the context of a smooth 
crossover regime at high and low temperatures. Such a study is of particular interest since extensive 
studies in lattice QCD over the past few years have led to the consensus that leans firmly in favor of a 
smooth crossover quark-hadron transition in the early universe. We showed that such a quark-hadron 
transition occurs and results in decreasing the effective temperature of the quark-gluon plasma during 
the process of the quark-hadron transition. The generic behavior of the temperature of the early 
universe in such a scenario is similar to that of a first order phase transition, although the differences 
in the energy should be taken into account. 
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